Finite-width effects in the near-threshold ZZZ and ZWW production at ILC 
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We calculate the cross-section of the near-threshold off-shell ZZZ and ZW + W~ production at 
the International Linear Collider taking into account their instability and the principal part of NLO 
corrections. The calculations are performed in the framework of the model of unstable particles 
with smeared mass-shell. We show that the contribution of the finite Z/W and H widths (their 
instability) is large in the Higgs resonance range and should be taken into account in the Higgs 
boson searches at future colliders. 



A great amount of work has been done so far in pre- 
cision tests of the Standard model (SM) including mea- 
surements of gauge boson, top quark masses and widths 
at LEP II [l|, l2| and Tevatron [3j , and very recently at 
much higher energies at LHC (see e.g. Ref. [1]). Due to 
clean environment and energies well above the electro- 
weak (EW) scale, future linear colliders would provide 
important tools for high-precision investigation ofgauge 
bosons and Higgs physics in the SM and beyond [5|, [6| . 

The multiple production of the gauge bosons is crucial 
for probing gauge boson (and Higgs) self-couplings, and 
thus for testing the non-Abelian structure and EW sym- 
metry breaking of the SM. The processes of two- (ZZ and 
W + W~) and three-boson (ZZZ and ZW+W-) produc- 
tion are of major importance as they give a direct infor- 
mation on trilinear and quartic vector boson couplings. 

Triple couplings of the neutral (Z and 7) and charged 
(W ± ) EW bosons, which were measured at LEP II [JQ 
and Tevatron Q, demonstrated a good agreement with 
the SM prediction within a few percent [l[ . For this pur- 
pose, the NLO EW factorizable corrections and finite- 
width effects (FWE) in the off-shell boson pair produc- 
tion i.e. e + e~ — » V*V* — > 4/ are very important, 
especially, in the near-threshold energy region (see e.g. 
Refs. [10, Lll I ) . However, corresponding higher-order cal- 
culations in the framework of traditional perturbation 
theory (PT) are rather cumbersome for 2 — » 4 processes 
as require evaluation of a few thousands one-loop dia- 
grams, and various schemes for automated loop calcula- 
tions are practically applied [H, [l| . 

Triple massive gauge boson (ZZZ and ZW + W~) pro- 
duction processes can be utilized to probe quartic gauge 
couplings and anomalous couplings in the Higgsstrahlung 
process (see, e.g. Ref. jLll - llq ]). At the moment, these 
processes being intensively studied in literature [17H20| . 
Typical leading-order contributions are shown in Fig. [T] 
In this work, we are primarily concentrated on the three- 
boson production processes in the Standard Model at 
ILC as the simplest case. 

Generally, FWE in a multiple gauge bosons produc- 
tion are closely connected to their instability, so they 
are usually referred to as the unstable particles (UP). 



FIG. 1: Diagrams for e + e" -» ZZZ, ZW + W~ processes. 



Near-threshold production of the unstable particles, as 
a rule, is accompanied by large FWE, which must be 
taken into account in analysis of corresponding observ- 
ables 0, 



In addition to the standard PT approach 
to FWE analysis, based on the stable particles approxi- 
mation (SPA), where UP instability is accounted for by 
higher-order corrections, various approximation schemes 
are practically applied in the literature, namely, semi- 
analytical approximation 2l[ 22 1, improved Born approx- 
imation [2^] , asymptotical expansions of the cross section 
in powers of coupling constant [24[ , fermion loop scheme, 
etc (see, also Ref. [25|, [26j and references therein). All 
above mentioned methods are based on the traditional 
quantum field theory of unstable particles [2l|. At the 
same time, there are some alternative approaches for UP 
description such as the effective theory of UP [27| , mod- 
ified perturbation theory [281 and the smeared-mass un- 
stable particles model [23, |30( . 

The main feature of the FWE is the "smearing" 
(fuzzing) of the threshold. In the standard treatment, 
this effect is described by taking into account all virtual 
states of UP, i.e. by its off-shellness. So, the cross sec- 
tion er(e + e~ — > VV) is defined as the cross section of 
inclusive four-fermion production cr(e + e~ — > 4/) in the 
double-pole approximation 31], which selects only dia- 
grams with two nearly resonant V bosons and the num- 
ber of contributing graphs is considerably reduced. Such 
an approximate description is usually realized with the 
help of the dressed UP propagators. 

In order to describe FWE in triple boson produc- 
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tion we have to consider full cascade process e + e — > 
ZZZ, ZW + W~ -> J2f 6 /i where the instability of the 
off-shell bosons is described by Breit-Wigner propaga- 
tors. So far, full NLO calculations were performed in the 
case of on-shell bosons only in the stable-particle approx- 
imation, i.e. without taking into account FWE and only 
for light SM Higgs boson with masses Mh = 120 GcV 
and 150 GeV [18l42dl |. The corresponding calculations in 
the traditional perturbation theory are rather cumber- 
some, as requires the complete set of a few thousands of 
one-loop diagrams. 

In the off-shell case, one encounters a very complicated 
problem. Exact NLO 2—^6 matrix elements would re- 
quire analytical evaluation of many tens of thousands 
loop diagrams, and are not available at the moment, but 
they are very important for boson FWE and Higgs con- 
tribution studies in vicinity of the threshold and Higgs 
resonance with the biggest cross section. 

In this paper, we describe FWEs in the 
near-threshold triple boson production e + e~ — > 
Z*Z*Z*, Z*W+*W-* -> E/ 6 / within the framework 
of the smeared-mass unstable particles (SMUP) model 
developed in Ref. [29|, [30| . In Ref. [32] the conception 
of the mass smearing as the main element of the 
SMUP model was successfully tested by comparison 
of its predictions with LEP II data on the boson-pair 
production e+e" -> Z*Z*,W+*W-* -> £/ 4/ total 
cross sections and Monte-Carlo simulations including full 
next-to-leading order matrix elements. In the framework 
of this model, the off-shell vector bosons Z, are 
treated as unstable particles, and the smearing of their 
mass shell effectively accounts for all-order propagator 
type corrections; the principal part of other factorisable 
corrections can effectively be taken into account as 
suggested in Ref. [32| , 

In this work, we present the total (inclusive) cross- 
sections of the off-shell ZZZ, ZW + W~ bosons produc- 
tion including FWE and the principal part of NLO (fac- 
torizable) corrections coming from the initial state radi- 
ation (ISR) and fussy-mass-shell unstable vector bosons 
in the framework of SMUP model. We found a large 
finite- width effect for Higgs masses Mjj > 2Mw,2Mz, 
when the contribution of Higgs diagrams becomes dom- 
inant and very sensitive to the Higgs boson width. So, 
this effect should be taken into consideration and can be 
applied as an auxiliary tool in Higgs boson searches at 
future colliders. 



BASICS OF THE SMUP MODEL 

In the off-shell ZZZ and ZW + W~ boson production 
in order to take into account FWE in standard way we 
need to consider full cascade process with six-particle 
final state and intermediate-state unstable Z and W 
bosons (see Fig. [2]). 



There is another way to tackle the issue. The SMUP 
model [2^, [3(| provides the possibility to treat Z and W 
bosons as final state particles and simultaneously to take 
into account their instability and obtain the finite-width 
effects correctly by smearing their mass shells [U l33j |. 
In this Section, we give a short description of the prin- 
cipal elements of the SMUP model and its advantages 
compared to the standard treatment (for more detail re- 
view of the SMUP model, see Ref. [3(1 EH and references 
therein) . 

The UP wave function in the framework of SMUP 
model is given by 



$a(x) = / $ a (x,p,)uj((i)d(i, 



(1) 



where <& a (x, A 1 ) is the standard spectral component which 
defines a particle with a fixed mass squared to 2 = fx 
in the stable (fixed- mass) particle approximation (SPA). 
The weight function w(/z) is then accounts for the self- 
energy interactions of the UP with vacuum fluctuations 
and decay products. This function includes all the infor- 
mation about UP decay properties (its instability) and 
describes the smeared ("fuzzed") mass-shell of the UP. 
The "fuzzing" of the UP mass shell is then caused, on 
the one hand, by quantum-mechanical instability accord- 
ing to the time-energy uncertainty relation and, on the 
other hand, by stochastic interactions of the UP with the 
electro- weak vacuum fluctuations. 

Then, the (anti) commutative relations for the UP field 
operators have an additional (5-function in the "smeared" 
UP mass (33 

$t(q,l/)]± = 8(ji-i/)6(k-q)6 a f>, (2) 

Here, subscripts "±" correspond to the fcrmion and bo- 
son fields. The presence of S(fi— p!) in Eq. @ means that 
the acts of creation and annihilation of the unstable par- 
ticles with different p do not interfere. So the quantity /1 
has the status of the physically distinguishable value of 
the UP mass squared to 2 . 

In the model under consideration, the transition am- 
plitude of the UP decay $ — > <px<p2 directly follows from 
Eqs. ([T]) and (J2|), and can be written as [30( 



A(k,p)=u J (,i)A st (k,p), 



(3) 



where A st (k,p) is the corresponding amplitude in the 
SPA, which is calculated in the standard way to a given 
order of the Perturbation Theory. From Eq. §3§ it follows 
that the differential (in UP mass squared fi) probability 
of the transition is dP(k,p) = p(fi)\ A(k, ii)\ 2 dp,, where 
p(/j>) = \uj(p,)\ 2 is probability density of mass parameter 
p = to 2 . This function is induced by multiple UP inter- 
actions with collective (in our case, EW) vacuum fluctua- 
tions of self-energy type and with decay products (vertex- 
type corrections) . In general, it is of non-perturbative na- 
ture and can be modeled in various ways. In this work, 
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we use the Lorentz distribution function 13 C 



Pip) 



Imn(fi) 



\}i - m 2 (p)] 2 + [Imll(^)] i 



(4) 



where II(fc 2 ) is the conventional vacuum polarisation 
function, m 2 (p) = + Rell(/i) with bare UP mass mo- 
The distribution ^ is of the Breit-Wigner type and ac- 
counts for the electro-weak quantum fluctuations of UP 
mass shell analogous to ones leading to a dressing up of 
the full propagator of an off-shell particle in the conven- 
tional quantum field theory. It naturally appears in the 
effective theory of UP as the most suitable one in the 
high-energy processes [3(| , and was tested before against 
LEP II data on off-shell ZZ and W + W~ near-threshold 
production [3^ |. 
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FIG. 2: Schematic illustration of the exact factorization of 
the total ee —> VVV —> 6/ cross-section in the SMUP model. 

As a general feature, such a model allows to consider 
the UP as the final state particles and their instability is 
then included by a convolution of decay widths or cross 
sections with UP mass probability density p(m 2 ) in each 
UP leg due to exact factorisation property which was 
proven previously in Ref. [34| . As a consequence of UP 
mass-smearing effect, this factorisation drastically sim- 
plifies the calculations. Now, we wish to apply the for- 
malism of SMUP model for prediction of the finite-width 
effects in the triple off-shell (ZZZ and ZW + W~) bosons 
production. 



MASS-SHELL SMEARING EFFECTS IN THE 
TRIPLE BOSONS PRODUCTION 

The processes e+e~ -> ZZZ; ZW + W~ at the tree 
level (leading order) in SPA are described by the set 
of diagrams represented in Fig. [T] Here, ZZZ produc- 
tion is described by nine diagrams with topologies a, &, 
whereas ZW + W~ production is given by sixteen dia- 
grams a, e. The first subset of diagrams a has a res- 
onant character at Mjj > 2Mn/,2Afz, and plays signif- 
icant role in the Higgs and FWE contributions. Com- 
plete NLO corrections to on-shell bosons production are 
described by additional few th ousands diagrams, which 
were calculated in Refs. 
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Due to exact factorisation property, the SMUP model 
3(1 34 1 allows to represent the total cross-section of the 



the factorized triple-convolution form (as schematically 
illustrated in Fig. HJ 



dm 1 dm 2 dm 3 a(s;m 1 ,m 2 ,rn 3 ) (5) 
xp(m\)p{m 2 2 )p(ml), 



where a(s; m 2 , m 2 , m 2 ) is the Born-level cross-section as 
a function of different bosons masses squared m 2 , i — 
1,2,3 in the SPA, p(m 2 ) is the probability distribution 
of the boson mass squared in ith leg given by Eq. 2] 
The Born ZZZ and ZW + W~ production cross-sections 
were calculated using FeynCalc v6.1 [13J as functions of 
smeared masses of Z and W bosons mf , and then their 
convolutions with p(m 2 ) over variable m 2 are performed 
numerically. In order to treat the poles in the boson prop- 
agators, which arise in the Higgs resonance region in the 
integration over the phase space of the boson pairs, we 
introduce the g 2 -dependent decay width of Higgs boson 
in the propagator |35| 



r H (q) = r%( q ) 



(?) 



1 H 



(<?), 



(6) 



where T s ^(q) is standard width of Higgs at the pole 



32j and other two terms account for the boson- 



pair channels of the Higgs decay. It should be noted here, 
that the exact Higgs width as a function of the momen- 
tum transfer scale (JS]) is very important since it deter- 
mines to a large extent the ZZZ and ZW + W~ produc- 
tion cross sections in the resonance regions, giving rise 
to the possibility of probing Higgs decay properties. If 
Higgs boson is heavier than 160 GeV and behaves as pre- 
dicted by the Standard Model, it decays predominantly 
into gauge-boson pairs and subsequently into four light 
fermions. And triple boson production at ILC can be 
rather sensitive to extra (anomalous) contributions to the 
HVV coupling from the new physics [16J. 

The approach under discussion has a close analogy 
with the convolution method [36| and the semi-analytical 
approximation 21, [22I . However, the status of these 



approaches are different 3^, 32, 33 1. In the framework 
of the SMUP model, the expression ([5]) directly follows 
from the UP smearing-mass conception, and the function 
p(m 2 ) describes the probability distribution in UP mass 
squared 2t| 3(|. Moreover, the definition of the unsta- 
ble particle field function ([1]) determines the strategy of 
taking into account the major part of the higher-order 
corrections in the near-threshold energy domain 32 , 33| ■ 



CROSS SECTIONS OF ZZZ AND ZW+W~ 
PRODUCTION AT ILC 

In the framework of SMUP model, the FWE were pre- 



viously studied in the boson-pair production e 



+ P - 



ZZ, W+W, ZH, Z-f at LEP II in Refs. [32j, 33]. In the 



inclusive process e + e — > ZZZ, ZW + W — > J2f®f m boson-pair production in vicinity of the threshold, the 
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NLO EW corrections are dominated by factorizable cor- 
rections to EW couplings, propagator-type (self-energy) 
corrections and the soft/hard initial and final state ra- 
diation while boxes, pentagons, etc can become relevant 
only at energies far from the threshold ones [26j |. 

As the next natural step, we apply the same strategy 
in analysis of the triple boson production. For this pur- 
pose, we perform the following consistency check - we 
compared our e + e~ — > ZZZ and ZW + W~ production 
cross sections in the SPA including only initial state radi- 
ation and NLO vertex corrections (renormalisation of the 
boson couplings) with full NLO results from Refs. 
in the 50 GeV range nearby to the threshold and got very 
close results within a percentage accuracy. This basically 
proves that non-factorisable, box and pentagon diagrams 
do not significantly contribute to the near-threshold pro- 
duction cross section, and can be omitted for our pur- 
poses. This is the only approximation we adopt in our 
calculations. 

Then we recalculated the on-shell ZZZ and ZW+W~ 
production cross-sections at fixed Mh = 120 GeV in SPA 
taking into account Higgs boson FWE. Again, the results 
coincide with ones reported in Ref. [II El. The contri- 
bution of Higgs diagrams and gauge boson FWE into the 
total cross-section at Mh = 120 GeV turns out to be rel- 
atively small. It is therefore very instructive to look at 
these contributions for heavier Higgs boson. 

The cross-section of the process e + e~ — > ZW + W~ as 
function of yfs at fixed Mh = 175 GeV (just above the 
latest Tevatron exclusion limit (33] ) is shown in Fig. |3] 
The Born cross sections in SPA and with taking into ac- 
count the Z/W FWE (due to the instability of gauge 
bosons in the final state) in the framework of SMUP 
model are represented by solid and dashed lines, respec- 
tively. The corresponding cross-section including the 
initial state radiation (ISR) corrections is given by the 
dashed-dotted line. These corrections together with the 
UP propagator-type corrections effectively taken into ac- 
count by UP mass-smearing effects in the framework 
of SMUP model are the main part of NLO corrections 
at considered energies [j^ 1 . One can see from Fig. [3] 
that the contribution of Z/W FWE is occurred to be 
quite large for relatively heavy Higgs boson in the near- 
threshold region. 

In Fig. |H the Born cross-section is shown as function 
of Higgs mass Mh at various fixed energies. Solid lines 
represent the total cross-section while the dashed lines - 
the part given by Higgs-less diagrams only. This figure 
illustrates rather strong dominance of the Higgs contribu- 
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FIG. 3: Cross-section of the process e + e~ — > ZW + W~ at 
fixed Higgs mass Mh = 175 GeV, given at leading order in 
SPA (solid line), including gauge bosons FWE (dashed line) 
and with taking into account both FWE and ISR corrections 
(dash-dotted line). 



tion due to Higgs resonance in the W + W~ phase space. 
It is worth to note here that the gauge boson FWE and 
smearing-mass effects naturally decrease at larger ener- 
gies, and get practically negligible at \fs > 450 GeV. 
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1 The major part of the vertex EW corrections at the threshold can 
be effectively included by taking coupling constants at the Mz 
mass scale, while contribution from the renormalisation group 
evolution is small at energies close to the threshold 1331 . 



FIG. 4: Cross-section of the process e + e~ — ¥ ZW + W~ as 
function of Mh at different energies. Contributions of the 
Higgs-less diagrams only are given by dashed lines. 

Analogously, the cross-section of the process e + e~ — > 
ZZZ as function of y/s at fixed Mh = 195 GeV close to 
the H — > ZZ threshold is given by Fig. [5l with the same 
notations as before. In Fig. [6] the Born cross-section is 
shown as function of Higgs mass at different fixed ener- 
gies. Again, similarly to the ZW + W~ case, the contri- 
butions of gauge bosons and Higgs FWE is large in the 
near-threshold region close to the Higgs resonance. 

In Conclusion, we would like to notice that using the 
argument that the mass-smearing conception realized 
in the smeared-mass unstable particles model is in the 
good agreement with the LEP II experimental data on 
the near-threshold boson-pair production, we applied the 
same ideas to off-shell ZZZ and ZW + W~ boson produc- 
tion at linear collides. The approach under consideration 
significantly simplifies the calculations with respect to 
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FIG. 5: Cross-section of the process e e — > at fixed 

Higgs mass Mh = 195 GeV. Notations here are the same as 
in Fig. E] 
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FIG. 6: Cross-section of the process e e — ► as function 

of Mh at different energies. 

the traditional one due to the exact factorization prop- 
erty. Explicit calculations of the triple-boson production 
cross sections demonstrate rather strong dependence of 
the gauge boson finite-width effects on the Higgs bo- 
son mass. In the near-threshold energy domain such 
effects are large and comparable with the initial state 
radiation corrections when the Higgs- resonant contribu- 
tion is significant, and small when the Higgs contribution 
is negligible. The Higgs-resonant contribution into the 
total cross-section is strongly dominant and have well- 
defined signature at Higgs masses above the Higgs decay 
threshold M H > 2M W for ZW + W~ production and at 
M H > 2M Z for ZZZ production. 
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